Universal depinning force fluctuations of an elastic line: 
Application to finite temperature behavior 
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The depinning of an elastic line in a random medium is studied via an extremal model. The latter 
gives access to the instantaneous depinning force for each successive conformation of the line. Based 
on conditional statistics the universal and non-universal parts of the depinning force distribution 
can be obtained. In particular the singular behavior close to a (macroscopic) critical threshold is 
obtained as a function of the roughness exponent of the front. We show moreover that the advance 
of the front is controlled by a very tenuous set of subcritical sites. Extension of the extremal model 
to a finite temperature is proposed, the scaling properties of which can be discussed based on the 
statistics of depinning force at zero temperature. 



I. INTRODUCTION 

The depinning of elastic interfaces allows to describe 
the phenomenology of a variety of physical phenomena 
such that the motion of a wetting contact line on a disor- 
dered substrate Q, the propagation of a fracture front 
in a heterogeneous material or the advance of a magnetic 
wall in a thin layer The depinning transition can be 
seen as a non equilibrium critical phenomenon (see Ref. 
[J HL El for recent reviews). The system is driven by 
an equivalent external force (the magnetic field for the 
domain wall, the stress intensity factor for the fracture 
front...) which plays the role of a control parameter. 
The richness of the phenomenology arises from the com- 
petition between the disordering effect of the quenched 
pinning potential and the ordering (smoothing) effect of 
the elastic forces acting on the front. The nature (short 
range/long range) of the latter will directly affect the 
universality class of the depinning. Moreover depending 
on the strength of the disorder, the motion front will ei- 
ther be collective (weak pinning) or consists of successive 
avalanches (strong pinning). In the following we focus on 
the latter situation. 

To study the pinning/depinning transition, a specific 
algorithm has been introduced under the generic name of 
"extremal model" 0,00; El- It consists in adapting the 
driving force so that only one site (the weakest one) can 
depin at a time. In so doing, the system is compelled to 
stay at the edge of the critical state (the front is driven at 
an infinitely low velocity). Thus the evolution consists in 
a series of equilibrium positions, because the driving force 
is adjusted to meet this condition. It thus can be seen as 
an ideal quasi-static driving[ll|. However, the price to 
pay in such a description, is that the dynamical aspects 
of the propagation are not accounted for (see Ref. 01 
for details on this subject). In a more complete picture 
either inertia or viscosity would have to be introduced. 
As dynamics is scarifies, time metric is not included in 
the "extremal model" description. What is preserved is 
simply the ordering of successive configurations. 

Close but below threshold, at finite temperature, an 
Arrhenius activation mechanism can allow for a steady 



"subcritical" propagation. Note that this situation differs 
from the creep regime at very low forcing which has been 
extensively studied in recent years |13lll4lll5lll6l |. In the 
latter, the driving force is very close to zero and thermal 
activation may allow back and forth motion along the 
direction of propagation. Close to threshold however, we 
may neglect the probability of a backward motion. More- 
over in the context of fracture or wetting, chemical reac- 
tions may produce an irreversible motion of the front. In 
the following we will not consider the possibility of front 
receding. The activated mechanism thus turns out to be 
relatively easy to implement as an extension of the ex- 
tremal model, where the latter is recovered in the limit 
of zero temperature. The motion consists in a succession 
of fast moves interrupted by long trapping events which 
require a thermal activation step to overcome an energy 
barrier. If the temperature is very low compared to the 
trap depths, the first site to depin will be the weakest 
one. In this limit of a vanishing temperature the system 
will thus naturally follow an extremal dynamics. The 
degeneracy of the time metric can however be clarified 
by studying the zero temperature limit. This descrip- 
tion is valid for a subcritical forcing, close to threshold. 
This forcing introduces a finite correlation length. How- 
ever, the thermal activation itself introduces another fi- 
nite correlation length. The competition between these 
two length scales will finally control the scaling property 
of the propagation. 

The aim of this paper is thus a careful study of the 
effect of temperature on both the dynamics and the crit- 
ical properties of the depinning fronts. A key parameter 
in the finite temperature depinning is the ratio between 
the thermal energy and the trap depths all over the front. 
The latter directly depend on the fluctuations of the de- 
pinning forces at zero temperatures. The external force 
needed to depin the front from a blocked configuration 
is actually a highly fluctuating quantity. The knowledge 
of the distribution of these depinning forces thus appears 
to be an essential ingredient in the study of depinning at 
finite temperature. 

The first part of this paper is thus devoted to the 
study of the depinning forces at zero temperature. We 
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show in particular that the distribution of dcpinning 
forces presents some universal features close to the crit- 
ical threshold. A criterion characterizing the transition 
between low and high temperature behavior requires ad- 
ditional information about the distribution of trap depths 
along the front. Again we show that this distribution 
presents some universal features. In the second part we 
introduce a Monte-Carlo like algorithm allowing to drive 
the front under constant forcing at finite temperature. 
We briefly present results concerning the two limiting 
cases (at low and high temperatures) and finally focus 
on the transition between the two regimes. The role of 
the two competing correlation lengths governed respec- 
tively by the driving force and the temperature will be 
discussed in this last section. 



II. DEPINNING FORCE DISTRIBUTIONS 
A. Description of the model 

At each time step, the front is represented as a sin- 
gle valued function h{x), with the additional assumption 
that its slope is small. The mean orientation of the front 
is along the x axis, while it propagates along a perpen- 
dicular direction y. For numerical simulations, the front 
is discretized on a regular grid of size Lwith periodic 
boundary conditions: hi = h(xi). 

Time is also discretized and incremented by one unit 
at each elementary move of the front. This time is thus 
a simple way to order the successive events, but it does 
not correspond to a physical time. Additional informa- 
tion has to be introduced to describe the off-equilibrium 
motion, and thus decide whether viscosity or inertia, or 
activation processes control the dynamics. The strength 
and weakness of the extremal approach is that only suc- 
cessive static positions are described, while the transi- 
tion between these conformations becomes part of the 
postulated rules of the model but do not stand for a real 
dynamic. Thus in this first part, t should not be misinter- 
preted as a real time. However, in the finite temperature 
section, we will come revisit this question and see how the 
extremal dynamics can be recovered with a more realistic 
dynamics. 

As introduced above, the motion of the front driven 
at an external force F depends only on the competition 
between a local trapping force j(x, h(x)) and an elastic 
interaction F el (x). We now specify them in more details. 

The distortion of the front due to the random envi- 
ronment induces elastic restoring forces. In the case of 
wetting the latter are the capillary forces. Using a small 
slope approximation, the elastic force contribution is lin- 
ear w.r.t. h(x), F el {x) = J G(x-x')h(x')dx'. Depending 
on the physical situation considered, these interactions 
can be short or long ranged. In the case of wetting the 
kernel G presents an algebraic decay G(r) oc r~ 2 up to 
the capillary length. In a Hclc-Shaw experiment how- 
ever, the cut-off scale of the capillary forces is given by 



the width of the cell and the Green function can be ap- 
proximated by G ~ 5" , the second derivative of the Dirac 
distribution <5. These short ranged interactions are then 
well described by a Laplacian term. In the following we 
focus on the latter situation with F el cx —d 2 h/dx 2 1 or in 
the discretized version, Ff = — + 2hi — 

The randomness of the environment can be introduced 
in two respects: first in the spatial distribution of the 
traps, second in the distribution of the trap depths. No 
correlation is considered here and these two quantities en- 
tirely characterized by their statistical distribution. Note 
that the same critical behavior is obtained as soon as ei- 
ther the trap positions or the trap depths are random. 
Either one of these two quantities can be a constant with- 
out changing the universality class of the model, i.e. all 
critical properties remain unaffected. We have performed 
a number of different numerical simulations changing the 
shape of the two distributions, and the only changes 
which are observed concerns small scales in space and 
time. In the following we will consider that the trap 
depth j(x,h(x,tj) and the front advance at the active 
site are uniformly sampled between and 1. 

In order to study the dcpinning force fluctuations, we 
drive the system with an extremal dynamics. This con- 
sists in adapting the external force so that only one site 
can depin at each iteration step. For a given external 
force F, the criterion for depinning at a particular posi- 
tion i on the front is 



F>F i (t)= 7i (h i (t))+Ff({h j (t)}) 



(1) 



One can thus naturally define the extremal site (i*, hi*) 
such that F c (t) = min^F^t) = 73* + F?j is the minimal 
external force to apply so that at least one site of the 
front can depin. The dcpinning criterion for a particular 
conformation of the front is thus 

F > F c (t) = min [ 7i (M*)) + F? ({hj (t)})} (2) 

The front depinning force F c (t) is in our simple version 
totally controlled by the front geometry. The extremal 
dynamics simply consists in tuning at each iteration step 
t the external force at exactly the depinning force of the 
current front conformation : F(t) = F c (t). From the dy- 
namical point of view it corresponds to drive the front at 
vanishing velocity with a very stiff spring. The system 
thus remains constantly at the edge of the critical state. 
The front explores a large number of different configura- 
tions to which correspond fluctuations of the depinning 
force F c (t). By definition the critical threshold F* is the 
minimum constant force to apply so that the front can 
advance indefinitely, it thus corresponds to the maximum 
of the front depinning force F c over all conformations: 

F* = maxF c (i) = max [min [n{hi{t)) + F? l ({hj(t)})] 

From the numerical point of view, the extremal dynam- 
ics finally amounts to iterate the following sequence:^ 
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FIG. 1: Standard deviation a(Ax) of the height difference be- 
tween points separated by a distance Ax for different system 
sizes L = 2 8 ,2 9 ,2 10 . The rescaling x -> x/L, a -» cr/L c_1 
with £ = 1.25 allows to capture the finite size effects. 



identify the extremal site i*; ii) advance the front at po- 
sition i* by a random increment; in) update the trap 
depth, on the site i* corresponding corresponding to 
the new front position h(i*); iv) update the elastic cou- 
plings Ff on the front to account for this local advance. 
The most time consuming step is the identification of the 
extremal site. Numerically efficient implementation can 
be used to reduce the computation cost of each iteration 
to log 2 (L) operations because of the short range nature 
of the interaction (note that in case of long range interac- 
tions, introducing an ultrametric distance along the front 
|l7j allows to reach the same numerical efficiency without 
changing the universality class of the model) . 



B. Scaling properties 

Let us summarize the expected scaling behaviors for 
an elastic line driven at a constant force F close to (but 
below) threshold F* . 

• Below threshold F* the length of the typical ad- 
vance £j_ scales as £j_ oc (F* — F)~ v± - 

• The correlation length along the front scales as £m oc 
\p* _ p\-v\\ _ 

• The relaxation time t of a line segment of length £ 
scales as r oc £?; z is called the dynamic exponent. 

• At threshold, in the steady state, the front con- 
formation is self-affine, characterized by a rough- 
ness exponent £. This can be shown by study- 
ing the average power spectrum of the front which 
scales as /c -1-2 ^. Similarly the wavelet transform 
also reveals the same exponent 01- I n rea l space, 
the typical width w (r.m.s. height) of the front 
over an interval Ax scales as w(Ax) oc (Ax)^, 
when £ is smaller or equal to 1. This is observed 



for some long-range elastic kernel. In our case, 
where the elastic kernel is the local curvature, we 
will see below that the £ exponent is larger than 
unity. In this "super-roughening" case, the scal- 
ing is anomalous [19| and the previous relation has 
to be corrected to w(Ax) oc L^~ 1 Ax. On Fig. 
we show numerical evidence for the validity of this 
"super-rough" scaling with the value £ = 1.25 for 
the roughness exponent. 

• Close to the critical threshold the motion consists 
of successive localized avalanches of lateral size Ax. 
The size distribution of these jumps follows a power 
law up to the correlation length £|| . Below £ii , the 
probability of observing jumps of size Ax > I scales 
as P(Ax > £) oc t~ x ~ A . In the framework of an ex- 
tremal dynamics driving, the driving force is not 
constant but is adjusted at the current depinning 
force F c (t). An avalanche triggered at a force F 
thus corresponds to a sequence of depinning events 
such that the current depinning force remains be- 
low the driving force F c (t) < F. The distributions 
of avalanches can thus be directly derived from the 
fluctuations of this depinning force F c (t). A di- 
rectly related quantity is the probability P(Ax\t') 
that after a sequence of r' events the depinning site 
has moved a distance Ax. We have 



P{Ax\t') oc Ax" a <; 



Ax 
Zap 



(4) 



where 4>{u) oc u a for u <C 1 and <f>{u) « cste for 
u 3> 1. In the framework of extremal dynamics, 
z' corresponds to a dynamic exponent (see Ref. 
p"2l l2(il ] for details on the relations between z' and 
the genuine dynamic exponent z introduced above. 
In the following we restrict the study to extremal 
dynamics and we use the notation z for the sake of 
simplicity. 

Simple scaling relationships immediately derive from 
these properties and allow to reduce the number of inde- 
pendent exponents. 

• Below threshold the front advance is confined to a 
region of length £11, and in the case of a roughness 
exponent £ > 1 the mean advance over this region 



is of order £j 
leads to 



oc 



£|ii c_1 oc (F- F*y v WLi- 



This 



V \ = V\\ = V 



with an unusual size-dependent prefactor. Note 
that this property is only true for a roughness ex- 
ponent £ > 1, otherwise u± = fii£. 

• Over a region of extent £|| , the typical macroscopic 
curvature is of order £jj -2 , and the latter scaling 
should balance the fluctuation of depinning force 
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FIG. 2: Distribution of distances between active sites corre- 
sponding to depinning events separated by a time delay be- 
tween At = 1 and At2 5 . After rescaling all distributions fall 
on a master curve. This rescaling has been obtained with 
the value 2 = 2 for the dynamic exponent. The dashed line 
corresponds to a power law of exponent a — —2.95 



observed in order to depin such a segment: AF oc 



£, Therefore^ 



= 2 



(5) 



• Let us recall that time only counts individual mo- 
tions, and is not physical. Therefore the dynamic 
exponent is an improper denomination. Neverthe- 
less, in this section, we keep this denomination 
since it refers to its common usage in statistical 
models. 

The mean advance of a segment of length £ scales as 
L c_1 £ and the size of an avalanche of lateral exten- 
sion £ writes r oc ^ 2 L^ _1 which leads immediately 
to 



z = 2 



(6) 



whenever ( > 1. Note however that even if the 
dynamic exponent z — 2 appears to be super- 
universal as soon as £ > 1, the relation between 
t and £ becomes system size dependent. Oth- 
erwise, for ( < 1, the same argument leads to 
z = 1 + CE3, E^- 
All exponents but a which characterizes the avalanche 
behavior can then be directly derived from the roughness 
exponent £. An simple scaling relationship can however 
be established in the context of the avalanches pi). Let 
us start from the probability P(Ax\t) that after a se- 
quence of r events the depinning site has moved from a 
distance Ax. The distance Ax is nothing but the sum 
of all successive jumps Ayi occurring at steps i 6 [1,t]. 
The distribution of these jumps being P(Ay) oc (Ay) _a . 
Assuming no time correlation we can apply a generalized 
central limit theorem for the sum variable Ax. The value 
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FIG. 3: Sequence of depinning forces F c (t). The dashed 
line corresponds to the value of the critical threshold F* = 
rnaxt F c {t) 



of the exponent a being slightly below 3, the limit dis- 
tribution is not Gaussian but a Levy distribution C a -i 
that exhibits at infinity the same power law behavior as 
the original power law distribution:£ a _i(Aa;) oc (Ax)~ a : 



P(Ax,t) 



-£a- 



(T) 



which leads to the scaling relationship a = 1 + z. In the 
present case of Laplacian elasticity, z = 2 and the pre- 
diction a = 3 slightly overestimates the value obtained 
by numerical simulations a « 2.95 (see Fig. This 
good agreement can be interpreted as the quasi absence 
of temporal correlations. Note however that in the case 
of long range elasticity where the role of temporal corre- 
lations is expected to be higher, this scaling relationship 
becomes a — 2 + ^ which is far above the measured values 
a ~ a where a is the exponent of the elastic redistribu- 
tion function. In the latter case, temporal correlations 
become more significant. 



C. Scaling and universality of depinning forces 
fluctuations 

On Fig. [3] we display a sequence of depinning forces 
F c (t) observed over 1000 steps of an extremal dynamics 
simulation and the corresponding distribution P(F C ) in 
Fig. 0] The upper force value of this distribution corre- 
sponds to the critical threshold F* . In the following we 
will argue that close to the critical threshold the distri- 
bution of depinning forces exhibits a universal behavior: 



P(F C ) cx (F* - F c y 



(8) 



Let us come back now to the series of depinning events. 
The depinning force F c fluctuates and can be described 
by a statistical distribution. The maximum value of the 
latter gives the value of the critical threshold. In the 
framework of a simulation this threshold F* is obtained 
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FIG. 4: Distribution of depinning forces, p(F c ). 

as the maximum depinning force over the ensemble of de- 
pinning events. Restricting ourselves to a finite system of 
lateral extension L, our determination will be biased by 
finite size effects. This invites us to introduce a simple 
prescription to get a very accurate determination of F*. 
Conditioning the distribution of depinning force by the 
distance d between consecutive active sites, we can built 
distributions P(F c ;d). We expect these distributions to 
converge to a Dirac distribution centered at F* as d di- 
verges since we probe only critically pinned configura- 
tions. This effect is directly observed on Fig. where we 
superimpose the original distribution of depinning forces 
P{F C ) and the contributions corresponding to increasing 
jump sizes d. We observe that, except for small jumps the 
successive distributions keep a similar shape up to scaling 
factors. In the following we check that these conditional 
distributions can indeed be rescaled onto a unique master 
curve using the scaling: 

i / F* — F \ 

nF,d) = - b ^{- F ^) (9 ) 

As a direct consequence of such a scaling we should in 
particular observe that the first and second cumulant of 
these distributions i.e the difference F* — (F c ) and the 
standard deviation both scale as d~ b : 

6F c (d) = F* — (F c )(d) occT 6 , (10) 
a Fc (d) = V(F?)(d)-(F c )2(d)ocd- b (11) 

This scaling directly leads to a linear dependence between 
(F c (d)) and <jF c (d). On Fig. El we check indeed that for 
large values of the jump size d, this linear behavior is 
obeyed. By extrapolation of this linear relationship for 
an infinite system for which the width ap^d) cancels, we 
obtain a precise determination of the critical threshold 
F*. We obtain in the present case 

F* = 0.392 ± 0.001 (12) 

Using the latter value we plot on Fig. 0the evolution 
of 8F c (d) and crp c (<i) versus d in logarithmic scale. The 



FIG. 5: Distribution of depinning forces conditioned by their 
distance to threshold, p(F c ,d), and the resulting global dis- 
tribution P(F C ) (bold). The larger the distance d between 
consecutive active sites, the narrower the distributions and 
the closer to the critical force. 
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FIG. 6: The linear relationship between the mean depinning 
force and its fluctuation allows to extrapolate the value of the 
depinning force for an infinite distance between consecutive 
active sites i.e. the depinning threshold F* for an infinite 
system. A numerical fit gives F* = 0.392 ± 0.002. 

observed linear behavior corresponds to power laws of ex- 
ponent b = 0.75. The value of the latter exponent can be 
understood since the depinning force F c directly derives 
from the local curvature. The latter can be estimated us- 
ing the elastic kernel and the roughness of the front. In 
the present case we have SF el (d) oc d c ~ 2 where C ~ 1-25 
is the roughness exponent of the depinning front. This 
estimation 6 = 2 — £ « 0. 75 is consistent with our numer- 
ical data. Moreover we can identify 1 jb to the exponent 
v characterizing the divergence of the correlation length 
close to threshold and we recover the scaling law of Eq. 
©: 

C + - = 2 (13) 
v 

Using this scaling we check on Fig. [SI that the succes- 
sive conditional distribution collapse onto a single uni- 
versal master curve. The knowledge of the distribution 



of distances between successive depinning sites allows us 
now to derive the behavior of the force distribution close 
to threshold. 



P(F*-F C ) = J P(F c ;e)p(£)de 



(14) 



(F* - F c )£- 



= (F* -F e )~ l+av ~" / u 1 - av+v - 1 i){u)du 



oc (F* — F c ) 



(a-l)iy-l 



We thus obtain a universal behavior of the distribution 
of depinning forces F c close to threshold as presented in 
Eq. E| with 



fi = (a — l)v — 1 ~ 1.6 



(15) 



Note in addition that when the distance d is small com- 
pared to the correlation length £, x = d/t; <C 1, the de- 
pendence of the P distribution w.r.t. (F* —F c ) should be 
independent of d, so that P(F C ; d) and P{F C ) share the 
same behavior (see FigJSJ . Hence, we expect ip( x ) ~ x ^ 
for x <C 1, in good agreement with the behavior shown 
on Figure H3 

At this point it is important to clarify that the distri- 
bution P by itself is not universal and depends on the 
local definition of the front advance and the trap depth. 
The details of the random distributions are absorbed in 
the small time and scale behavior. In contrast collective 
effects contribute to large jumps in the active sites and 
'J is universal as well as the singular behavior of P close 
to threshold. 
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FIG. 7: The mean distance from the depinning force to the 
depinning threshold and the depinning force fluctuations cal- 
culated for given distances between consecutive active sites 
obey a scaling law with the latter distance. The straight lines 
plotted as a guide to the eye correspond to power laws of 
exponent 0.76 and 0.77 consistent with the expected value 
2 - C » 0.75. 
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FIG. 8: After rescaling the distributions of depinning forces 
corresponding to distances d = 4, 8, 16, 32 collapse onto a 
unique universal master curve. 



D. Subcritical sites along the depinning front 

In Fig. [5] we show the distribution of all individual 
depinning forces Fi = ji — Ff along the front obtained 
from extremal dynamics. We see clearly that two dif- 
ferent populations of sites can be distinguished. Above 
threshold F*, Q(F > F*) follows some kind of trun- 
cated Gaussian distribution. We also note that there is a 
small but significant number of sites which carry a force 
smaller than threshold, F < F* . Moreover, the decay of 
frequency is very abrupt right below F*. We will argue 
that this decay can be rationalized as a power-law of the 
difference (F* - F c ), Q(F < F*) oc A(L)(F* - F)~ e . In 
the insert of Fig. [5]we plotted the proportion of subcriti- 
cal sites Q(F < F*) against the system size. This evolu- 
tion appears to follow a power law Q(F < F*) cx L D_1 
with D w 0.35. We show below that D can be regarded 
as a fractal dimension of the set of subcritical sites along 
the front. 

In the Appendix, we will introduce a mean-field ver- 
sion of the model where such a distribution can be com- 
puted analytically. In this simple variant of the depin- 
ning model, the elastic force is not transfered from the 
active site to its nearest neighbors but to two sites chosen 
randomly along the front. Although a quantitative com- 
parison can not be established because of the mean-field 
nature of the solvable model, on a qualitative ground, the 
main features of the force distribution are recovered. On 
FigEB one can distinguish again between a Gaussian-like 
distribution for depinning forces above the critical thresh- 
old and a singular part below the threshold. Moreover it 
can be shown that the amplitude of this singular part is 
system size dependent. 

The "subcritical" sites which carry a force smaller than 
threshold, play a very important role, as we now try to 
clarify. At each time step, the extremal site corresponds 
to the minimum force and always remains smaller than 
F* . Thus the active site is always part of the subcritical 
sites. Its depinning induces a local unloading and an 
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FIG. 9: Distribution of individual depinning forces Fi = 7i — 
Fi along the front. The subcritical part of the distribution 
(F < F*) is size dependent. In the insert we see that the 
relative weight of this subcritical part follow a power law of 
exponent 0.65 with the size L. This suggests that the set of 
subcritical site is fractal with a dimension D = 1 — 0.65 = 0.35 



additional re-loading on its neighbors. We note that a 
subcritical site has a chance of being activated wherever 
it is located, while a site carrying a force larger than 
threshold can only become active if one of its neighbors 
become active itself. These sites can thus be regarded as 
a population of "precursors" of depinning. 

In order to determine the distribution of the depinning 
forces along the front we use the same tools as above. We 
first show that the distribution of F{ along these sites 
depends on their distance to the extremal one as a power 
law: Q(Fi < F*) oc \x t — Xi*\~ a . We observe then that 
the conditional distribution Q(Fi < F*; \xi — Xi*\ = d) 
obeys the same scaling as the one studied above for the 
distribution of depinning forces. This allows finally to 
estimate the power law behavior of the distribution of 
the depinning forces of this population of subcritical sites 
close to threshold. 

On Fig. 1111 we plot in logarithmic scale the mean dif- 
ference to threshold (F* — i^; \xi — Xi*\ = d)) for sites 
at a distance d from the extremal one i*. We obtain 
the same scaling as for the distribution of the depinning 
force F c conditioned by the distance between successive 
depinning sites: 



{F* - Ff, \xi —Xi*\ = d)) oc d 



2-C (16) 



and again we see on Fig. Elthat we can rescale all condi- 
tioned distributions onto a single master curve provided 
the distance d is large enough. 



P(F*-F i = f>0; 



= d) = d-v[fd- (17) 



On Fig. E| we show that the subcritical sites are not 
homogeneously distributed along the front. We plot in 
logarithmic scale the distribution of distances d along the 
front between the extremal site and any other subcriti- 
cal site. The straight line obtained in logarithmic scale 



indicates that the latter obeys a power law: 
n (\xi* —Xi\ = d;F < F*) oc dT c , 



(18) 



with c w 0.68, what indicates that the subcritical sites 
form a fractal set of dimension D = 1 — c w 0.32, which 
is consistent with the value D « 0.35 obtained above for 
the evolution of the population of subcritical sites against 
the system size. 

The above scaling relationships can now be used to 
estimate the distribution of depinning forces along the 
front close to threshold: 



Q(F* -F i = f>0) 

A{L) f Q(F* - Ff,£)n(£)d£ 



(19) 



= A(L) J l~ (j) \ {F* — Fi)* 

= A(L) (F* - Fi)~ 1+cv ~ v 
oc A(L) (F* - Fi)- 5 



r c de 



with e = 1 + vD « 1.45. The size dependent prefactor 
can be estimated using the fact that for a front of size L 
the typical cut-off of the elastic force is L~ x l v . Rewriting 
the total number of sub-critical sites along the front then 
gives: 



Af(L) = LA(L) 



[F* — Fi)~ 



LA{L)L-^-^ 
A(L)L 1+D 



(20) 

(21) 
(22) 



The fractal dimension of this population of subcritical 
sites being D, we have in addition N(L) oc L D which 




FIG. 10: Comparison between the result of a direct numerical 
simulation (symbols) and the previous analytic expressions 
(continuous curves). The system size was set to L = 2 s , 2 10 
and 2 12 . As above we see that only the subcritical part is size 
dependent 
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FIG. 11: The difference 8F (current force to threshold) of 
precursory sites located at a distance d from the extremal site 
scales as a power law of d. The straight line corresponds to a 
power law of exponent 0.75 «2-(. 
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FIG. 12: After rescaling the distributions of subcritical depin- 
ning forces Fi conditioned by the distance d to the extremal 
site i* , collapse onto a single master curve. 



Q(F* - Fi 



f > 0) cx i (F* - Fi)' 



1 



D 



2-C 
(23) 

We see thus that the depinning is controlled by the evo- 
lution of a very tiny fraction of subcritical sites. These 
sites are distributed fractally along the front and can be 
regarded as precursors of the avalanches to come. Close 
to threshold the depinning force associated to these sites 
is power law distributed with an exponent depending 
both on the roughness exponent of the front £ and the 
fractal dimension D characterizing the population of sub- 
critical sites. 



III. FINITE TEMPERATURE DYNAMICS 

We now consider the case of depinning at finite tem- 
perature. The transition loses its critical character. As a 
matter of fact, even at arbitrary low external driving, the 
front keeps moving. The velocity is however extremely 
low, most of the time being spent to thermally depin 
from blocked conformations. This creep motion has been 
intensively studied in recent years from the experimen- 
tal H and theoretical points of view 

El El El El- it 

has been argued that at very low forcing, the velocity fol- 
lows a stretched exponential behavior against the driving 
force: 



v(F, T) cx exp 



F* 



T 
To 



(24) 



where m = (2£ — l)/(2 — (). To our knowledge, the 
anomalous force scaling predicted in Eq. (|24|l has not 
been reproduced so far by numerical simulations |22l |23| . 

Note again that such a behavior is only expected to 
happen at very low forcing. This means in particular 
that the macroscopic motion of the front results from a 
balance between microscopic backward and forward mo- 
tions. 

In this section we focus our study on a slightly differ- 
ent situation, namely, a driving force in the vicinity of the 
critical threshold. This allows to neglect the probability 
of a backward motion. At finite temperature the mo- 
tion consists of a series of rapid finite moves interrupted 
by long pinned stages from which the front can only es- 
cape by thermal activation. Two limit cases can actually 
be distinguished regarding the thermal depinning of a 
blocked conformation of the front. At very low temper- 
ature, the front most likely depins at the extremal site. 
Thus the sequence of depinning events is not affected by 
the introduction of a temperature. The computation of 
the time spent by the front to escape any blocked con- 
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FIG. 13: Distribution n(d) of distances d of subcritical sites 
to the extremal site. The straight line indicates a power law 
of exponent 0.32. The data being obtained with a logarithmic 
sampling of d, this corresponds to a power law n(d) oc d~ c , 
c^l- 0.32 » 0.68 
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FIG. 14: Distribution of subcritical depinning forces F < 
F* along the front. Close to threshold the distribution of 
subcritical sites follow a power law of exponent —1.75. 



in this quasi extremal dynamics. The knowledge of the 
depinning force distributions P(F C ) thus serves to quan- 
tify the velocity of the front. The second limit case corre- 
sponds to high temperatures. During a transient, all sites 
are free to advance; then the curvature of the front be- 
comes comparable to the temperature and the dynamic 
of the front is rather well described by an annealed dy- 
namics. This induces a scaling of the front roughness 
different from the one obtained at zero temperature: we 
observe a roughness exponent £y w 0.5 instead of the 
zero temperature value £ ~ 1.25. 

In the following we first present the simple Monte Carlo 
algorithm used to run the dynamics of the front at finite 
temperature. The latter is obtained by allowing acti- 
vated depinning along the front. We then give results 
obtained in the low and high temperature limits. We 
finally present a criterion allowing to discriminate be- 
tween these two regimes. It appears actually that beside 
the "static" correlation length £ of the critical transition 
another characteristic length of thermal origin £x plays 
a key role in the problem: below £t the front follows 
a T = dynamics but above £t the activated depin- 
ning become dominant and the front follows an annealed 
Edwards- Wilkinson like dynamics [2^. 



A simple Monte-Carlo algorithm for depinning 
at finite T 



Let us come back to the definition of the dynamics. 
Consider a point i along the front. At zero temperature 
the criterion for depinning at this particular site is: 



Di = F e 



■Ft - 7i > 



(25) 



Note here that for the sake of simplicity we only con- 
sider the fluctuating part of the random traps, the con- 
stant part, call it V being implicitly absorbed in the ex- 
ternal driving force: F ext = — T. The values of 
this effective driving force may thus be arbitrarily low 
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FIG. 15: Dependence of the mean waiting time on the tem- 
perature after rescaling by the dominant Arrhenius factor. 
The straight line indicates a power law behavior of exponent 
2.3 to compare with the value (a — l)u « 2.6 obtained in Eq. 
1321 for the low temperature dynamics. 



and even negative without violation of the hypothesis of 
microscopic forward motion mentioned above. 

At finite temperature we can estimate the local depin- 
ning time ti via an Arrhenius term: 



1 



• exp 



n 



l 

Ma 



1; exp 



-Di 
T 



(26) 

where Vq corresponds to a natural frequency of vibration 
at the microscopic scale. 

From the numerical point of view we can use the fol- 
lowing algorithm. At time t, all points satisfying the cri- 
terion of propagation l|25[) are advanced simultaneously 
up to the next traps. The waiting times Tj corresponding 
to the points j that do not satisfy the criterion are com- 
puted from their respective distributions. Points such 
that the waiting time is below the time unit are also ad- 
vanced. The elastic forces are then updated and time is 
incremented. If no point along the front has been ad- 
vanced, the location of the first site to thermally depin 
and the associated waiting time are drawn from the Ar- 
rhenius distributions. Again the selected point moves 
forward, the elastic forces are updated along the front 
but the time is incremented by the waiting time spent 
in this blocked conformation. Let us now be more ex- 
plicit about the determination of the waiting time and 
the location of the depinning site in the case of a blocked 
conformation. 

The probability Pi(t)dt that the site i is the first to 
thermally depin in the time interval [i, t + dt] can be 
written: 



m) = 



i 

— exp 

Ti 
1 

— exp 







r^ exp ( 









(28) 



The probability r, that the activated depinning event 



takes place at this particular site i is thus independent 
of time 



Pi(t)dt 



(29) 



In a similar manner, the probability P ac t(t) that the first 
activated depinning event takes place at time t is 



1 ft 
Pact(t) = — exp I - — 



(30) 



Note that the time t and location i appear as indepen- 
dent variables. When the front is pinned, the algorithm 
thus consists of drawing a waiting time from the distribu- 
tion P act (i) and to choose the site of depinning according 
to the weights rj . The characteristic times involved in Ti 
depend on the distribution P(Fi) of individual depinning 
forces along the front. In the following we focus on the 
two limiting cases of low and high temperature. 

At low temperature, when the front is pinned the wait- 
ing time is dominated by the characteristic time of the 
extremal site. The sequence of depinning events is iden- 
tical to the one of the extremal dynamics. The time 
evolution directly results from the knowledge of the dis- 
tribution of depinning forces P{F C ) of the extremal. At 
high temperature, the probability of depinning becomes 
comparable for all sites along the front. 



Low Temperature behavior: 
dynamics 
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FIG. 16: Mean wavelet coefficients obtained for the depin- 
ning front driven at Fext = 0, (F* — F ex t = 0.569) for various 
temperature. The straight lines indicate the two limit be- 
haviors of roughness exponents £ = 1.25 and (t = 0.5 corre- 
sponding to the critical state for temperatures close to zero a 
"SOS-like" behavior for higher forces. System size N = 4096, 
32 10 6 iterations. 



observe a simple Arrhenius behavior for the mean depin- 
ning time with a power law correction in the ratio tem- 
perature/distance to threshold. On Figure fl~5l we show 
the scaling of this waiting time. After rescaling by the 
Arrhenius factor we observe indeed a power law depen- 
dence on the temperature. Numerically we find for the 
exponent v(a — 1) = 2.3 to compare with the expected 
value v(a — 1) = 2.6 with £ = 1.25 and a — 2.95. 



We consider here the case of very low temperature. 
The loading external force F ext is constant below the 
threshold F*. At a given conformation of the front, all 
sites satisfying the criterion of propagation are advanced, 
the elastic forces are then updated. If the front is blocked 
the only site to depin is the extremal one. The time r 
associated to an iteration is drawn from an Arrhenius 
distribution 



p(r) 



1 



■ exp 



1 



max 



1; exp 



F,, 



T 



(31) 

The knowledge of the distribution of depinning forces 
P(F* — F c ) gives us directly the distribution of waiting 
times in blocked configurations)^. The mean waiting 
time (t) to depin from a blocked conformation for an 
external constant loading F ext is 



(t) = / di^ c exp 

J Fext 



l] " F, r ' t '\A{F*-F c Y{2,2) 



exp 



T ) 
F* — F t 

r r ext 
f 



(33) 



where we used the condition T <C F* — F ext to trun- 
cate the integral in the numerator. In this regime we 



C. High Temperature behavior: 
Edwards- Wilkinson like dynamics 

At high temperature, all growth probabilities tend to 
be equal, which generates an independent evolution at 
each site. The result is simply a white noise morphol- 
ogy, with ever increasing slopes and curvatures. This 
regime can only be a transient one. At long times (and at 
lower temperatures) , a bias between the growth probabil- 
ity and the curvature will be felt and hence, this limit will 
correspond to an annealed model known as the Edwards- 
Wilkinson (EW) model[24[> which indeed give rise to self- 
affine fronts with a roughness exponent ( = 0.5. We thus 
expect to see a cross-over from ( = 1.25 at small tem- 
perature, to £ = 0.5 at larger temperature. This is pre- 
cisely what is observed in the numerical simulations as 
shown in Figure IT^l which represent the scaling behavior 
of the front roughness for a same driving force but differ- 
ent temperatures. Here we measured along the front the 
averaged wavelet coefficients oj(a). For a self-affine front 
of roughness exponent £, these coefficients are expected 
to scale as cj(a) oc a 1+2 ^ with the length scale q[l8|. In 
log — log scale we observe a scaling behavior with an ap- 
parent exponent decreasing from £ w 1.25 toward £ w 0.5 
for growing temperatures. 
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FIG. 17: Time evolution ol the interface width at high tem- 
perature (T = 0.1) with a low forcing F* — F ext = 0.39. 
The time and height coordinates are rescaled by L Z EW and 
jJ,ew respectively. The value zew ~ 2 and (ew = 0.5. The 
evolution of the width is well described by the power law 
w(t) oc t 0EW with (3ew = 0.5 (dashed lines). These criti- 
cal exponents are consistent with a growth of the Edwards- 
Wilkinson type. 



To check that this high-temperature behavior indeed 
belongs to the EW universality class, we show on Fig. El 
the time evolution of the interface width after rescaling 
of the time by L z and IS with z = 2 and £ = 0.5. As 
expected in a EW growth model we observe a power law 
behavior with an exponent = 0.25. 



D. Transition from low to high temperature: 
thermal characteristic length 

Let us now consider intermediate temperatures. We 
can actually establish a quantitative criterion to discrim- 
inate between low and high temperature behavior. 

As stated above the probability of depinning of a 
pinned configuration derives from the distribution of de- 
pinning forces along the front. These depinning forces 
can be ranked from the smallest (the extremal site) to 
the largest one. The extremal site is denoted i* and the 
associated depinning force is F c . The second smallest is 
ii and the corresponding force is Fi 2 . The typical depin- 
ning time t* is: 



1 = E 1 = -( 1 +E- 

T* ^ TV TV. ^ TV, 



(34) 



3=2 "I 



where 



tv» ( F ext - F c F ext - Fa 
-=exp f — — j=<*l> 



Fjj ~ F c 
T 



(35) 

We see here that the relative weights corresponding to 
the different sites do not depend on the level of external 
driving but only on the difference of force to the extremal 



FIG. 18: Difference between the two smallest depinning 
forces F2 — F c along the front against the difference F* — F c . 
The data are well described by power law of exponent e = 0.7. 



one. It is thus clear that as soon as F c — F^ ^S> T the 
waiting time r* is dominated by the extremal time Tj» . 
The ranking allows to restrict to the criterion: 



F c - Fi, > T 



(36) 



The previous section gave us the knowledge on the dis- 
tribution of depinning forces along the front, namely: 



Q(F* - Ft) <x ^(F* - Ft] 



(37) 



We may thus estimate the typical difference between the 
two smallest depinning forces along the front: 



*2 1 



F*-F c 



Yp{F* — F)~ e dF = I ^4 

(38) 

where 1/L D in the integral is a normalizing factor of the 
subcritical force distribution and L D at the right hand 
sides is the number of subcritical site. A Taylor expan- 
sion around F c leads to 



(F c - F l2 ) w (F* - F c f 



(39) 



We see in Fig. El that this scaling is nicely observed 
with the value e = 1.7 consistent with the estimation ob- 
tained above for the distribution of subcritical depinning 
forces. The extremal value of the depinning force F c is 
a fluctuating quantity during the course of the propa- 
gation. In the first part of this paper we saw that the 
difference F* — F c can be associated to a typical dis- 
tance d between consecutive depinning sites by the scal- 
ing d oc (F* — F c ) ~ v . The criterion for extremal depin- 
ning during the propagation can thus be characterized 
by a distance £t below which the activated site is the 
extremal one and above which the activated site is ran- 
domly chosen along the front: 



£t oc T 



(40) 
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FIG. 19: After rescaling the abscissa by the characteristic 
length £t <x T~ u ^ with v = 1.33 and e = 1.4, and the ordi- 
nate by £y with £ = 1.25, the mean wavelet coefficients ob- 
tained for the depinning front driven at the constant driving 
force F em t = 0, (F* — F ex t = 0.569) for different temperatures 
collapse onto a single master curve. 



In Fig. we presented the scaling of the front rough- 
ness obtained for a constant forcing at various temper- 
atures. As expected, the scaling evolves from the zero 
temperature behavior £ = 1.25 to an annealed behavior 
C = 0.5: the various curves are distributed between these 
two limit cases. However, rescaling the lateral length by 
the temperature dependent characteristic length £t, we 
can collapse all results onto a single master curve (see 
Fig. I19fl where the small arguments correspond to the 
critical behavior and the large arguments to the annealed 
behavior. 



IV. CONCLUSION 



among this set. The latter can thus be regarded as pre- 
cursors of the future avalanches. All of these sites are 
characterized by a sub-critical depinning force, namely 
the force F needed to make them move is below the crit- 
ical threshold F* . Again the distribution of the depin- 
ning forces associated to these precursory sites presents 
a universal character: P(F) oc L~ 1 (F* — F)~ e , again 
composed of a universal form when conditioned by the 
distance to the active site. 

In a second part we extended the study of the de- 
pinning transition at finite temperature. The front was 
driven at a constant below threshold driving force. A 
simple Monte-Carlo algorithm was developed to simulate 
the activated dynamics. 

Two limit cases were considered. At very low temper- 
ature or at small scale, the extremal site is by far the 
most probable site to depin and all critical properties 
are preserved. In this limit, the finite temperature then 
mimics the extremal dynamics. Assuming an Arrhcnius 
behavior for the depinning from any individual trap we 
are able to associate a physical time to this quasi ex- 
tremal dynamics. The knowledge of the depinning force 
distribution P{F C ) oc (F* — F c )^ then allows a quantita- 
tive estimate of the front velocity. We recover a simple 
Arrhenius dependence with a temperature power law de- 
pendence. At high temperature or at large scales, all 
sites can equally be activated. This annealed dynam- 
ics has been shown to be well described by an Edwards- 
Wilkinson growth model. The transition between these 
two regimes can be estimated using the ratio of the tem- 
perature to the gap between the two smaller depinning 
forces. We showed that this gap can be naturally asso- 
ciated to a temperature dependent characteristic length 
along the front £t oc T^ 1 ^" '. Below the characteristic 
length £ the front adopts naturally a critical dynamics, 
at larger scales the annealed dynamics dominates. 



We studied the depinning of an elastic line in strong 
pinning conditions. 

In a first part we focused on the critical behavior at 
zero temperature, based on an extremal dynamics to per- 
form numerical simulations. Beyond the scaling proper- 
ties of the front roughness and the avalanche dynam- 
ics observed at the depinning transition our main results 
concern the statistics of the depinning forces. At each 
conformation of the front we can associate the minimum 
force F c (t) needed to depin the front. The distribution 
of these depinning forces present a singular behavior in 
the vicinity of its maximum value, the critical threshold 
F* is P{F C ) ex (F* - F C ) M where the critical exponent 
fi directly depends on the roughness exponent £ of the 
depinning front. Moreover this distribution relative to 
large jumps d in the activity adopts a universal form for 
all values of F. 

Moreover it appeared that the critical dynamics con- 
centrates on a very tenuous set of sites (termed subcrit- 
ical) along the front, in fact a fractal structure. Be- 
yond neighboring effects, the critical site is always chosen 



APPENDIX A: A SOLVABLE MEAN FIELD 
DEPINNING MODEL 

We present in this appendix a mean field variant of the 
extremal depinning model studied in section II. We first 
note that we can set all thresholds to 0, and provided that 
the moves are still random, the behavior remains similar. 
In this context the extremal site i* simply minimizes the 
local elastic force: 

F c = Ffi = mmFf . (Al) 

i 

Once the extremal site has been determined the front 
advances up to the next trap. However we introduce 
here a change in the definition of the elastic forces. In 
the original Laplacian version Ff l = — + hi — 
the extremal site is unloaded proportionally to the local 
displacement just covered and conversely its two nearest 
neighbors are equally loaded. In the mean field version 
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we present here, the two sites loaded to keep the elastic 
force balance are chosen randomly along the front, hence 
the term mean field. 

Let us call x the gap between the local elastic force 
and the critical threshold, x = F el — /*, p(x) its dis- 
tribution density, and P{x) the cumulative probability 
P{x) = J p{y)dy (from which p(x) = —P'(x) results). 

For a system of size L, the smallest x has a density 
p s {x) and cumulative P s (x). We have P s (x) = P(x) L , 
and thus p s (x) — LP{x) L ~ l p{x). 

Ignoring spatial correlations, we can write a mean-field 
master equation for p{x) . 



dp(x) 
dt 



= -Ps(x) - 2p{x) 

+ (1/2) / p s (x-y)dy 



(A2) 



p(x + y)dy 

= -p a (x) - 2p(x) 

+ (l/2)(P s (x - 2) - P s (x)) 
+2(P{x)-P(x + l)) 

The terms of the right hand side correspond to the 
probability of occurrence of the following events i) x is 
the extremal value of the elastic force along the site; ii) 
Before the update, one of the two random sites to be 
loaded carried a force x; in) After update the advancing 
site carries a force x; iv) After update one of the two 
random sites carries a force x. 

In the steady state, we have 

-2P'(x)-P' s (x) = (l/2)(P s (x-2)-P s (x))+2(P(x)-P(x4 

(A3) 

In order to integrate this equation, we distinguish in 
the following between three different intervals. The re- 
gion 1 (x < 0) corresponds to the subcritical sites. The 
region 3 (x > 1) corresponds to deeply pinned sites which 
have no chance to become neither subcritical nor critical 
at the next depinning event. The intermediate region 2 
(x > 1) exchanges with the two other ones. For x < 0, 



we expect P(x) to be equal to 1 but a small correction 
vanishing as L diverges. P(x) = 1 — ^(x). We look for 
y L (x) = (l/i)$(x). Thus P s (x) = (1 - (l/L)<f>(x)) L -> 
exp(— $(x)). The equations corresponding to the three 
regions write: 

• Part 1: 

- P' s (x) = (1/2)(1 - P s (x)) + 2(1 - P(x + 1)) , (x < 0) 

(A4) 

• Part 2: 

- P'(x) = P(x) - P(x + 1) + (1/4) , (0 < x < 1) (A5) 

• Part 3: 

-2P'(x) = (l/2)P s (x-2) + 2P(x) , (x>l) (A6) 
and the solutions are respectively 



• Part 1: (x < 0) 

P.{x) 
P(x) 



~x z - 2x , 
{ 1 2(x + l) 
L x(x + 2) 



Part 2 and 3: (x > 0) 



P(x) = (l/4)(x-2) 2 



(A7) 
(A8) 



(A9) 



These analytical expressions are compared with nu- 
T))merical results on Fig. ^] Beyond the excellent agree- 
ment between analytics and numerics, the similarity be- 
tween Fig. El and ED shows that this simple mean field 
model allow to recover the main feature of the original 
model. The distribution of local depinning forces is com- 
posed of a size independent supercritical part and of a 
size dependent subcritical part which presents in addi- 
tion a singular behavior close to the forces approach the 
critical threshold. 
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